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Abstract: We explore the phenomenological predictions of a supersymmetric standard 
model, with a large extra dimension and unifying gauge couplings. The modified five 
dimensional renormalisation group equations make it possible to obtain light, maximally 
mixed stops, with a low scale of supersymmetry breaking and a low unification scale. 
This allows the fine-tuning to be lowered right down to the barrier coming directly from 
experimental lower limits on the stop masses. We also show that modifying the SUSY 
breaking pattern to obtain lighter stops at the high scale does not result in fine-tuning 
relaxation, and only RGE effects turn out to be effective in generating a lower fine-tuning. 
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1 Introduction 

The discovery of the Higgs boson, of mass rrih ~ 125 GeV, at the first LHC run [1, 2] 
and various null results in searches for super-particles (sparticles) appear to imply that the 
sparticles of minimal theories of supersymmetry (such as the MSSM) are likely to be out 
of reach of the LHC altogether. Further considerations regarding electroweak symmetry 
breaking and the naturalness aesthetic bring further doubt that perhaps supersymmetry 
is not a symmetry of nature after all. 

At this time it is also worth to consider non minimal models that share much of 
the well grounded theoretical elegance of the (four dimensional) MSSM: supersymmetry, 
gauge coupling unification, anomaly cancellation, a minimal matter content to achieve 
these, unification of matter representations, an explanation of the big hierarchy problem, 
various dark matter candidates. The model we consider is five dimensional and brings with 
it additional interesting features: the possibility to observe Kaluza-Klein states of gauge 
(and other) fields [3], the possibility to achieve the observed Higgs mass with sparticles 
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within reach of the next LHC run, and a much lower unification scale and supersymmetry 
breaking scale than is normally possible in four dimensions. 

In this paper we explore electroweak symmetry breaking in a particular class of five 
dimensional supersymmetric theories [4, 5]. We also wish to address the possibility to 
construct a ‘natural theory’ where stops are lighter than their first and second generation 
counterparts, and which is not spoiled by the renormalisation group effects that would 
usually make the 3rd generation similarly heavy at a low scale, after a long period of 
RG-running. Studying models that differ in how matter fields are located among branes 
and bulk opens the possibility to explore the effect different SUSY breaking patterns and 
modified five dimensional RGEs have on fine-tuning wrt a four dimensional theory. Un¬ 
fortunately modifying the breaking pattern to obtain more natural soft terms at the high 
scale does not give the expected fine-tuning relaxation. However we show that modified 
RGEs enable us to obtain light maximally mixed stops. This allows us to lower fine-tuning 
down to the barrier coming directly from observational lower limits on the stop mass. 

In our analysis we used renormalisation group equations outlined in [5, 6] and adapted 
a G-|—|- based spectrum generator originally intended for the (four dimensional) MSSM [7]. 
A similar modification may be carried out with any publicly available spectrum generator 
[8-10] The RGEs used in this paper may be found in [5] and further conventions in 
[6] and [15-19]. Eor earlier phenomenological studies of five dimensional theories see for 
example [20]. 

The outline of the paper is as follows: section 2 we outline the two base models that 
we wish to explore. In section 3 we explore the behaviour of the running parameters of the 
theories. In section 5 we look at naturalness and electroweak symmetry breaking as we vary 
the radius of the extra dimension. We also look at benchmark models of supersymmetry 
breaking in section 4. In section 4.3 we outline and explore a “natural” model in which the 
3rd generation are spatially located on a different brane to the first two generations and 
in which supersymmetry is gauge mediated directly to the first two generations but only 
indirectly to the 3rd, which ideally would allow for light stops. In section 6 we conclude. In 
appendix A we outline the implementation of spectrum generator and RG solver in C-|—|- 
code. 

2 The 5D-SSM+(F±) Model 

The first model that we wish to explore is a five dimensional supersymmetric theory with 
the field content outlined in table 1 and is pictured in figure 1 (left). In this model the 
Higgs fields gauge fields and additionally are bulk fields [21]. This matter 

content is necessary for the gauge couplings unification, as we shall explore further later. 
All five dimensional bulk matter fields are supersymmetric Hypermultiplets which due to 
even and odd boundary conditions lead to a four dimensional Ghiral multiplet as a zero 
mode of the Kaluza-Klein expansion: such details are well documented, for instance in 

^To date, five dimensional theories are one such class of models that cannot yet be explored using SARAH 
[11-14] although it can still be a powerful tool to determine the RGEs of the low energy four dimensional 
effective theory that the five dimensional theory runs to [5]. 
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Figure 1. Pictorials to represent the location of matter in the five dimensional model. In Model 1 
(left), all generations of matter live on a brane. In Model 2 (right), only the 3rd generation live on 
a brane 
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Table 1. The matter content of model 1. All superfields of chiral fermions live on a brane and all 
Higgs-type superfields and gauge vector fields live in the bulk. The superscript / = 1,2,3 denotes 
the generations. Neutrino superfields may be included straightforwardly. 


[15-19]. The second model we wish to explore is outlined in table 2 and pictured in figure 
1 (right). In model 2 only the third generation is located on a brane and the first and 
second generation are in the bulk along with the Higgs multiplets and fields. 

The superpotential for both models is given by 

W=YuU Hi - Yd d djf Hi - Ye e + piHuHd + fiP-F+ . (2.1) 

It would be very worthwhile to consider the generation of the term jlF~F~^ in the su¬ 
perpotential, although for this paper we will not need to consider it, and postpone that 
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Table 2. The matter content of model 2. 


to later work. We will now explore the running parameters of these two theories as one 
changes the scale of the extra dimension. 

3 Running parameters 

It is particularly interesting to understand and compare the behaviour of the various run¬ 
ning parameters of these theories compared to the more usual four dimensional MSSM. 
The behaviour of the various parameters as a function of renormalisation scale for model 1 
is pictured in figure 2. Of particular note is that unification happens much earlier if the size 
of the extra dimension is large [22] , than the usual four dimensional case. One also finds 
that the top Yukawa reduces rather significantly and becomes of similar order to the other 
Yukawa couplings near the unification scale. In addition one finds that even for initially 
vanishing A-terms the At term may become multi-TeV in value at the electoweak scale, 
which is encouraging from the perspective of obtaining the observed 125 GeV Higgs mass. 
It is also the case (bottom left) that the gluino mass can become much hearvier than the 
other gauginos allowing for the theory to still have a light bino and wino whilst allowing 
for a gluino above current exclusions. 

The first model may be compared with model 2 similarly presented in figure 3 and in 
table 2. In these figures it is notable that that gauge couplings quite nearly unify but the 
gauge couplings rise rather than fall, after the KK modes start to take effect in the RGEs. 
The Yt still decreases in value, although now rather interestingly the At becomes so quickly 
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negative that it can quickly overcompensate the effect of the gluino soft mass, and for very 
large radius, the At running may even return on itself. Again the wino and bino soft terms 
can be much smaller than that of the gluino, even starting from the same initial value. 



Log,o(fJ/GeV) Log,[,(fj/GeV) 


Figure 2. Model 1 running of gauge coupling constants gi{fi) (top left panel), 3rd generation 
Yukawa couplings (top right panel), trilinear soft terms (bottom right panel) and gaugino soft terms 
(bottom left panel) with compactification scales 1/R ^ 10^ GeV, 10® GeVdi 10^^ GeV, as a function 
of Logio(p,/GeY). In this example all soft terms were set to Msusy = ITeV at the unification 
scale (defined by gi = ( 72 ), except the trilinear soft terms [At) which were set to 0. 
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Figure 3. Model 2 running of gauge coupling constants gi{fi) (top left panel), 3rd generation 
Yukawa couplings (top right panel), trilinear soft terms (bottom right panel) and gaugino soft terms 
(bottom left panel) with compactification scales 1/R ^ 10^ GeV, 10® GeV© 10^^ GeV, as a function 
of Logio(p,/GeY). In this example all soft terms were set to Msusy = ITeV at the unification 
scale (defined by gi = ( 72 ), except the trilinear soft terms [Ai) which were set to 0. 


4 Supersymmetry breaking in benchmark models 

So far our exploration has been reasonably agnostic about how supersymmetry is broken, 
since the main feature of the models presented in the previous sections are their RGBs. In 
what follows we will simply refer to sets of RGBs we used, as models. 

There are however a number of ways that have been proposed for the parametrisation of 
supersymmetry breaking in a five dimensional scenario. In this section we wish to identify 
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these scenarios and look at their patterns of supersymmetry breaking which define their 
possible high scale spectra. 

4.1 Gravity mediation (CMSSM) 

Our first benchmark scenario is the simple CMSSM spectrum, however since easier genera¬ 
tion of A-terms during running is a key feature of five dimensional running, we will always 
take Ai = 0 case for which the difference between five and four dimensional theories is the 
most visible. This implies a very simple type of spectrum with just two free parameters 
Ml and mo: 

Mi = Ml , = m-o , = 0, (4.1) 

defined at the unification scale. 


4.2 Minimal Gauge Mediation (MGM) in five dimensions 

In gauge mediation there is an additional characteristic scale at which SUSY is broken, 
which for brevity we will labelled M. For the five dimensional RGEs to have an impact on 
the spectrum and to not simply be an effective four dimensional theory with a low SUSY 
breaking scale we wish that M is at least 0{1/R) and possibly nearer Munification- 

The soft terms in five dimensional GMSB, at the breaking scale, are then given by 



where F and M are the free parameters we will scan over. This paper is the first im¬ 
plementation of five dimensional GMSB soft masses [16, 17, 23-25], with five dimensional 
RGEs [5, 6]. In both model 1 and 2 we will take the supersymmetry breaking to be on the 
opposite brane to the matter, and both brane and bulk matter are essentially suppressed 
by the effect of the extra dimension, as in the above equation. 


4.3 Realising natural SUSY with GMSB in five dimensions (nMGM) 

The renormalisation group equations of model 1 may be used to explore a natural susy 
scenario as pictured in figure 4. In this model the 3rd generation is located on one brane 
and the 1st and 2nd generation on another, along with the supersymmetry breaking sector. 
The effects of supersymmetry breaking are mediated by gauge forces [26] (but one can also 
easily consider gravity mediation too in this context) and the result is that the 1st and 2nd 
generation and also the gauginos will receive normal (4D) GMSB soft mass contributions 
but the 3rd generation will be heavily suppressed [5, 6, 17, 25]. The soft mass matrix for 
squarks and sleptons takes the form 




/l 0 o\ 
0 10 
(^0 0 oy 


... 


(4.3) 


leading to an interesting natural SUSY spectrum of lighter 3rd generation squarks. This 
scenario suggests that natural SUSY softer terms are imprinted due to the ‘geometry’ of the 
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Table 3. Experimental exclusion limits used 

theory. We will consider such a natural spectrum in context of minimal gauge mediation, 
the resulting soft terms are similar to those in(4.2), however now only third generation 
sfermions are suppressed by In the text we will refer to this as an nMGM 

spectrum. Needless to say, a similar model may be constructed using brane to brane 
gravity mediation. It would also be interesting to discuss models with Hu and Hd localised 
alongside the 3rd families, however it would require a much more serious modification of 
the RGE’s of our Model 1 and 2, consequently we postpone that discussion to future work. 



Figure 4. Pictorial to represent the location of matter in the five dimensional model resulting in a 
natural SUSY breaking scenario with GMSB (nMGM spectrum) and model 1 RGBs. 


5 Electroweak symmetry breaking and naturalness 

One important feature of a model is whether its parameter space can accommodate elec¬ 
troweak symmetry breaking. Figure 5 shows regions in the parameter space of our models 
where the breaking does not occur or which violate direct detection bounds summarised 
in table 3 [27]. Exclusions corresponding to varying size of the extra dimension (including 
the AD case) are plotted together. For standard GMSSM and MGM boundary conditions 
Model 1 predicts rather standard spectra of sparticles quite similar to the AD case. How¬ 
ever Model 2 due to much lower gaugino masses compared to the M-terms allows us to 
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obtain very light stops and maximal mixing even despite A-terms vanishing at the unihca- 
tion scale. In fact for large R = 10“^ the peculiar shape of the CMSSM excluded region 
in model 2 comes from obtaining too light stops that would have already been observed. 

The MGM excluded region comes from the interplay between large scalar masses we 
obtain at the scale M when M = 1/R, and when they are generated during 5D modihed 
running between 1/R and M » 1/R. The minimal stop mass is obtained between these 
two situations and results in excluded part on the left hand side of middle row in Figure 5 
where the small stop soft mass fails to push to negative values and break electroweak 
symmetry. 

This is also visible in nMGM plot on the bottom row of Figure 5. However here the 
problem is more severe since is not suppressed by 1/{MR) at the SUSY breaking 
scale, and a much bigger part of the parameter space is excluded. For nMGM spectrum 
this problem appears also for very small l/(Mi?), because in this part of the parameter 
space the difference between Higgs and stop soft masses is the largest. These two effects 
lead to appearance of a window of allowed parameter space which is very interesting, since 
it is in that window, that we obtain the highest Higgs mass. 


5.1 Exploring naturalness in benchmark scenarios 


In MSSM-like theories, at a hnite loop order, electroweak symmetry breaking is radiatively 
induced. The up-Higgs soft mass is driven to negative values, leading to 


“ 2"^^ ~ + 1^1^ + ^ (tan/3 

At leading order, the running of this soft mass in four dimensions follows 

3 


~ (mq, + rriu^ + )Log 


Msusy \ 

Ms J 


(5.1) 


(5.2) 


In hve dimensional models the RGEs are rather different due to the power law contributions 
and one finds 


--+ '^Us + ^t) 


Log 


1 


RMs 


+ AIsusyR 


(5.3) 


One might have expected a significant contribution to fine tuning from the power law 
contribution. However four and five dimensional theories actually have similar fine tuning 
as the much faster power law contribution can dominate the running for only a very small 
range of scales if the spectra we are comparing are similar. And so the final amount of 
fine-tuning for a given scenario depends mostly on the resulting spectrum rather than 
on the amount of power law running. This is quantified in figure 6, where in numerical 
calculations we use a standard fine-tuning measure with respect to parameter a defined as 
follows^ [28-30] 


A 


a 


din Ml 
5 In a 


(5.4) 


2 


The numerical procedure used to calculate fine-tuning is detailed in appendix 


A.5. 
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Figure 5. Striped regions of the CMSSM (top row) MGM (middle row) and nMGM (bottom row) 
parameter space cannot accommodate electroweak symmetry breaking or are already excluded by 
direct searches. Left hand side plot shows results for model 1 and right hand side for model 2. Both 
show different sizes of the extra dimension and CMSSM shows the 4D case as well. 


Fine-tuning connected with a set of independent parameters Oj is then 
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(5.5) 
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Figure 6 shows resulting fine-tuning as a function of Higgs mass for different sizes of the 
extra dimension as well as the result one would obtain from 4Z) running, the top row shows 
results obtained assuming CMSSM-like soft terms (with Ai = 0), the middle row shows 
gauge mediated boundary conditions and the bottom plot shows the nMGM ones. 

The results in left panel show model 1 which gives a rather standard prediction despite 
power law contribution to running. However model 2 shown on the right hand side allows 
us to reduce fine-tuning very significantly. The reason are the gaugino masses that decrease 
during 5D part of the running (as shown in Figure 3). This protects the soft terms from 
the usual increase due to the heavy gluino. Since the A-terms do not grow proportionally 
to scalar masses we can easily achieve maximal mixing scenario for the light stops, and 
their direct detection bound is precisely what gives us the lower bound on hne-tuning we 
can see in model 2 with R = 10““^. 

The bottom plot shows nMGM result which turns out quite similar to MGM and 
CMSSM model 1 results. The reason for this is that in model 1 the least fine tuned 
results are those for which M » l/R. Thus the scalar masses are initially very small 
and have to be generated with modified running. Gonsequently the 3rd family part of 
the spectra are very similar. The correction introduced by nMGM relies only on larger 
subleading corrections to the Higgs mass from first two families and other Higgs sector 
scalars. Unfortunately hne-tuning price of these corrections is larger than their contribution 
to the Higgs mass and the results are slightly more hne tuned than those from standard 
MGM or CMSSM soft terms. 

A large qualitative difference between MGM and nMGM becomes visible for Higgs 
masses slightly higher than the observed one. This comes from the part of parameter space 
which predicts successful electroweak symmetry breaking in nMGM. As explained in the 
beginning of this section, the problem is a result of the exclusion appearing in nMGM 
for very small 1/MR. Where we cannot break electroweak symmetry because radiative 
correction to the unsuppressed soft Higgs mass coming from highly suppressed stop mass 
is to small, and the former never runs negative. This becomes visible for higher Higgs 
masses because very small 1/MR is the part of the parameter space where we obtain 
highest Higgs masses. Another very important feature of 5D models is the possibility to 
bring superpartner masses within the LHC reach for points predicting minimal fine-tuning. 
This is illustrated in Table 4 which shows spectra corresponding to lowest obtained hne 
tuning for nih = 125 GeV. 

6 Conclusions 

In this paper we explored the implementation of the hve dimensional renormalisation group 
equations of a number of supersymmetric extensions of the MSSM, into a full C-|--|- spec¬ 
trum generator, along with self energy corrections for the Higgs mass. 

Our key result is showing that modihed hve dimensional RGEs can result in spectra 
very different from the usual AD case. The is because in 5D the heavy gluino does not 
necessarily dominate running of other soft terms during power law running, as in our model 
2. Thus we can easily obtain maximal stop mixing and much less hne tuned spectra, even 
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Figure 6. Fine-tuning as a function of Higgs mass for different sizes of the extra dimension for 
models 1 (left hand side) and 2 (right hand side) with CMSSM (top row), MGM (middle row) and 
nMGM (bottom row) spectra as well as the AD results. 


with standard sets of soft terms at the SUSY breaking scale. This is also very interesting 
because in 5D models the least fine tuned spectra with correct Higgs mass can easily 
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Table 4. Masses of superpartners (in TeV) for spectra which minimize fine-tuning for mh = 
125 GeV 


predict soft superpartner masses within LHC reach, even for standard patterns of soft 
terms. Interestingly, this means the most interesting parts of the parameter space can be 
probed during next run of the LHC, which is not usually the case in 4D models. 

We explored models where the 1st and 2nd generation are in the bulk and a model 
in which the 1st and 2nd generation is on the same brane as the supersymmetry breaking 
sector and the 3rd generation is located on an opposite brane, resulting in a spectrum of 
stops lighter than other squarks. Obtaining lighter stop soft terms at the SUSY breaking 
scale did not result in a more natural spectrum. The reason is non negligible fine-tuning 
price of heavier first two generations and heavier Higgs sector which give only a subleading 
correction to the light Higgs mass. 

The final advantage is a low scale of unification of gauge couplings and a low supersym¬ 
metry breaking scale. And also much better unification of Yukawa couplings (especially in 
model 2) which gives hope for a very interesting five dimensional UV completion of such 
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models. 
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A Numerical procedure 

In this appendix we outline the implementation of the RG-solver and spectrum generator 
used in this paper. The numerical procedure we use is similar to the ones used in existing 
codes [8-10]. We work with quantities renormalized in DR and use renormalization group 
equations (RGE), to iteratively find low energy parameters for a given set of high energy 
soft terms. 


A.l Mz Scale 

At the scale Mz we include radiative corrections to couplings. We set Yukawa couplings 
using the tree-level relations 


rritV^ 


mf ,\/2 


11 Li V ^ 11 LhV ^ Tfl V ^ 

yt = —^ , Vb = -X , Vt = (A.l) 

V sm p V cos p V cos p 

where mt,mh,mr are fermion masses and v is the Higgs vev. At the first iteration we use 
physical masses and SM Higgs vev v ~ 246, 22 GeV. During subsequent iterations above 
quantities are renormalized in DR scheme and one-loop SUSY corrections are included. 
To calculate the top mass we use 2-loop QCD corrections [31] and 1-loop corrections from 
super-partners from the appendix of [32]. While calculating the bottom mass we follow 
Les Houches Accord [33], starting from running mass in MS scheme in SM Next 

applying the procedure described in [34] we find DR mass at Mz, from which we get MSSM 
value by including corrections described in appendix D of [32]. While calculating the tau 
mass we include only leading corrections approximated in [32]. We calculate the Higgs vev 
in the MSSM using 

5i + 3<7?/5 ’ 

where we include Z self interactions described in appendix D of [32]. To calculate gi , g 2 , 
g '3 in DR in the MSSM we use the procedure described in appendix G of [32]. 


-\/2 


A.2 RGE and scale 

After calculating coupling constants at the scale Mz we numerically solve RGEs [35], [36], 
to find their values at the scale M„, at which we include the soft breaking terms. Then we 
solve RGEs a gain to find soft terms, coup ling constants, tan f5 and Higgs vev v at the scale 
Mew SB = {^ewsb)bii^ {Mewsb)- At first iteration we take p = sgn(p) GeV and 

Bfj^ = 0 and run to the scale at which the above equation is fulfilled. 
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Figure 7. Schematic of the numerical algorithm. Subsequent 
steps are described in the appendix. 


A.3 Electroweak symmetry breaking 

In order to obtain correct electroweak symmetry breaking we use minimization conditions 
for the scalar potential to find new values of ^ and We include radiative corrections 
in these equations by the substitution 


mHu + — 

Vn, 




td 

ruHa H-• 

Vd 


(A.3) 


We include full one-loop corrections to tu and presented in appendix E of [32] and 
leading two-loop corrections [37-41]. Since these corrections depend on sparticle masses 
which in turn depend on the /i parameter that we aim to calculate, an iterative calculation 
is performed to obtain new values of /r and 

If the new values differ significantly from the ones obtained in previous repetition of 
the whole algorithm described above, we run back to the Mz scale and repeat the whole 
calculation once again. If however the values of /r and converged, we can move on to 
teh calculation of physical masses. 
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A.4 Calculation of physical masses 

To calculate physical masses we use only leading corrections described in [32] everywhere 
but the Higgs sector. In the Higgs masses calculation we use full one-loop corrections from 
[32] and leading two-loop corrections described in [37-41]. 

A.5 Fine-tuning 

After the calculation of the spectrum is finished, one has a whole set of parameters and 
couplings that predict correct electroweak symmetry breaking. In order to calculate fine- 
tuning we solve the RGBs from Mu scale down to Mew SB with one of the fundamental 
parameters a* changed slightly at the high scale Mu ■ Than at the scale Mew SB we 
recalculate the spectrum and use minimization conditions to calculate a new value of tan /3 
and to obtain our new prediction for which means that we calculate numerically the 
derivative in the definition of fine-tuning (5.4). We repeat that procedure for all parameters 
Ui and obtain our final result as a maximum of results obtained for each of those parameters 
(as in (5.5)). 
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